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I. INTRODUCTION 


The motivation for the work reported here came from two 
exploratory projects, a Rotary Detonative Wave Engine, 
(RDWE), described in Ref. 1 and a high speed Centrifugal 
Diffuser Test Device ‘(CDTD), described in Ref. 2. Both 
projects required the geometrical design of an axisymmetric 
generator for nearly-tangential transonic outflow. 

To improve gas turbine engine performance an increase 
in cycle pressure ratio is needed. One alternative to 
mechanical compression is to supply the compression 
thermally by detonative combustion. If the pressure rise 
produced in this manner is sufficient for the engine, no 
compressor and associated drive turbine are required and a 
simpler more efficient engine might result. 

As described by Monks [Ref. 1] , a proposed geometry 
for such an engine might include a rotor with outflow 
passages which are turned at the exit to produce nearly 
tangential flow. Figure 1 illustrates the proposed RDWE. 
Monks attempted to measure the torque on a static simulated 
rotor as a detonation wave propagated from the center 
outward. A reaction torque was produced from the nearly- 


tangential expulsion of gases from the rim. 
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The question then was, what would be the best design 
for the rotor passage to maximize the torque impulse? In an 
analysis included in Appendix A it was shown that tangential 
flow at the rim of the RDWE would yield the highest specific 
power and highest power conversion efficiency. A design was 
required therefore of nozzles which would direct the outflow 
almost tangentially. The nozzles would operate in an 
unsteady mode. 

The second project which motivated the present study 
involved the development of a test rig for radial 
diffusers. Such diffusers are used in centrifugal 
compressors to convert the kinetic energy of the flow from 
the rotor into static pressure. Erwin (Ref. 2] proposed and 
designed a so-called low speed Centrifugal Diffuser Test 
Device, ‘CDTD). The CDTD simulated the flow of air 
delivered to the diffuser from the compressor rotor using 
atatic swirling vanes. The facility would allow the 
evaluation of proposed new diffuser geometries and permit 
validation of computer analyses for diffusers operating in 
the fully subsonic flow regime. Figure 2 shows a schematic 
of the low speed CDTD. Vidos [{Ref. 3] presented results 
from the flow generator used in the CDTD and performed an 
analysis of the internal flow. As mentioned by Vidos, the 
design of centrifugal diffusers is presently based on 
experimental results for two dimensional and conical 


diffusers. Numerical methods are currently under 
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development for viscous, three dimensional, transonic flow 
Ca regime containing both supersonic and subsonic regions), 
with adverse static pressure gradients. Such methods are 
required to design optimum radial diffusers for transonic 
conditions. Equally, a high speed version of the CDTD was 
required in order to examine such diffuser flows under 
controlled conditions and to verify the emerging codes. 
Again, nearly tangential swirling outflow was required, this 
time in steady flow. 

The main interest here was therefore to design vane 
Shapes and the apparatus for a high speed CDTD. The design 
would closely parallel the general arrangement of the low 
speed CDTD. The high speed CDTD, had the unusual 
requirement of radial-outflow with convergency at the rim 
exit. The limitation on the design was primarily that the 
hardware could be machined easily and inexpensively in- 
house. Computer programs for the design of supersonic 
nozzles for radial-inflow turbomachines were found to exist, 
but most were highly involved and would need extensive 
machining processes to generate the blade shapes. 
Particular designs are discussed herein. A literature 
search indicated that computer programs for radial-outflow 
devices using simple vanes for two-dimensional isentropic 
flow of a perfect gas were not available. 

In the present report, the design of the blading for a 


high speed CDTD model and the general arrangement of the 
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apparatus are described. General methods and computer 
programs for the design of double circular arc and wedge arc 
blading for radial-outflow vanes are described in detail in 
appendices. The report also discusses a method of 
predicting the flowfield through a selected passage using an 


Euler code based on the Godunov method of solution. 
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II. APPARATUS DESIGN 


Factors considered in the design of the high speed CDTD 
fell into two categories; namely, limitations due to the 
available air supply system and piping and the requirement 
that the test section have dimensions appropriately scaled 


for the diffusers required to be tested. 


A. DESIGN CONSIDERATIONS 

The high speed CDTD was designed to operate as an 
intermittent blow-down wind tunnel. In order to use the air 
supply system available in building 230 at the Naval 
Postgraduate School, the test section flow area was dictated 
by the available mass of stored air and the minimum desired 
runtime of two minutes. Using proportions similar to those 
used in the low speed CDTD (see Figure 2), the radius of the 
swirl generator became four inches while the radius at the 
test section was five inches. The method used for 
calculating these dimensions is given in the following 
section. The difference between the swirl generator outer 
radius and the test section radius allowed for mixing of the 
wakes from the blading (or nozzles) used in the swirl 
generator. The awirl generator was limited to a length of 
two inches to meet the scaling requirements set by the test 


Section. Because of the larger pressure differential 
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required across the swirl generator in the high speed 
version, turbine like turning vanes were selected rather 
than the sheet metal nozzle construction used in the low 


speed device. 


B. APPARATUS DESCRIPTION 

Figure 3 illustrates the overall arrangement of the 
proposed CDTD design. The swirl vanes for the high speed 
CDTD are sandwiched between two circular plexiglas disks. 
Plexiglas disks were used to allow visualization of the flow 
through the vanes. The disks are mounted on an aluminum 
center pipe shown in Figure 3. The pipe provides double 
axial in-flow. Following the work performed by Feieresien 
C[Ref. 4] a transition section to turn the air from the axial 
to the radial direction was added in the form of a plate on 
a tie-rod on the pipe centerline. The design of the 
transition section is not complex if the area of the radial 
outlet is smaller than that of the total axial inlet. An 
acceleration through the transition section will then occur, 
leading to generally favorable pressure gradients. These 
consideration required a double-entry supply arrangement for 
the CDTD model. 

The piping is attached to the air supply system. Figure 
4 shows the configuration of the blades on the plexiglas 


disks and in relationship to the supply pipe. 
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C. DESIGN CALCULATIONS 

Vavra (Ref. 5] showed that in the design of a turbine 
cascade the efflux angle « can be calculated from the 
geometry of the cascade and that it is independent of the 
incidence angle i of the inlet velocity, see Figure 5. 
Vavra also noted that it is common practice in steam turbine 
design to use the relation 


© = arc cos ( a/ sa) C1) 


where “a” is the throat opening at the discharge and s is 
the blade spacing. Equation (1) applies when the trailing 
edges of the blades are of negligible thickness. [In 
practice the “wedge” angle at the trailing edge and trailing 
edge bluntness modifies this dependence. 

With the limitation of a two minute runtime and the 
available mass of stored air, Demo [Ref. 6] found that for 
the supersonic cascade rig in building 230 an area of 6.27 
square inches was needed at the nozzle throat where the flow 
would choke. This same area was taken for the total outlet 
throat area of the swirl vanes for the high speed CDTD. 

Referring to the following sketch, the area at the exit 
where the flow will choke, A*, is given by, 

Ax = L a 2Z@ (2) 


where L is the vane length, a is the throat width and Z is 


the number of vanes. 
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Figure II.1 Description of Conditions at Throat 


The blade spacing, #, is given by 
ase nds Zz ec D 
where D is the diameter at the vane exit. From equations 
(2) and (3), 
Ae = L (Ca/s) 1D. (4) 
In Ref. 7, the following conditions were calculated for 
a “full scale” high speed CDTD design to have a test section 
width of 1 inch at 25 inches diameter: 
M1271.5 at 831255°, at R 1212.5 inches 
V121392 ft/aec at To=520 °R, giving 
V 6421140 ft/sec and VRo21357 ft/sec and for 
a swirl vane length of 5 inches, Vpo5=182 ft/sec 
and Bo282.34°. An increase of swirl vane length 
to 5.713 inches gave VRo"160 ft/sec and 5283.30 
Using the above data for the "full acale” high apeed 
CDTD in equation (4), A#* = 38.5 inches2. Using a linear 


acaling to the model in which A* = 6.27 inches2, the scale 


factor is 
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(Aepsauy)* = (6.27/38.5)7% = .40. - (3) 


Using this sacale factor, the dimensions calculated for the 


high speed model are shown in the following sketch. 


= [O" OIA. 


{~—---——— _ |= 2" ——————_ | 


= ee 


fa on. 


Figure II.2 Cross Sectional View of the Test Section 


From equations (2), 
azz? Ae / L (6) 
therefore for a diameter of 8 inches we can predict a 
required throat dimension, and blade spacing using equation 
(3) for a selected number of blades. Results are given in 
table 1. The requirement then is to design suitable turning 


vanea to produce the required throat area and spacing. 
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III. BLADING DESIGN 


A. BACKGROUND 

In Ref. 5 Vavra gave a procedure to establish a blade 
profile for specified values of throat width a, spacing s, 
trailing edge thickness te, and inlet flow angle, « Figure 
6 illustrates these features. 

Vavra’s design was for axial turbine blades or vanes. A 
computer program was available to construct blading 
following Vavra’s methods (Ref. 5]. The method used 
circular arcs, and straight lines. Wedge angles and leading 
and trailing edge radii were input parameters to determine 
@ flow channel with suitable acceleration. No comparable 
method and program were found for radial outflow vanes. 

In 1968 Goldman and Scullin (Ref. 8] designed blades 
with inlet transition arcs, circular arcs and outlet 
transition arcs. However, these blades were not converging 
at the exit and were not designed for radial-outflow. 

Later Goldman and Vanco (Ref. 9] developed a computer 
program for the design of sharped-edged throat supersonic 
nozzles for isentropic flow. Their program was constructed 
for the radial-inflow case. The Goldman and Vanco design 
also involved the requirement of contour machining. 

There is a substantial amount of information on both 


radial-inflow and axial turbine blade design in the 
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literature C[Ref. 8,9,10,11, and 12 for example]. Yet none 
was found which dealt with the required case of radial- 


outflow. 


B. DESIGN METHOD 

The expansion process in a radial turbine differs 
appreciably from that in an axial turbine because of the 
radius changes. Vavra’s work (Ref. SJ] involved axial 
turbine design. The work presented here for radial-outflow 
vane design followed an appreciation of Vavra’s methods and 
possibly retraced steps of others who have previously sought 
the simplest geometries for machining purposes. 

The simplest shape, that of double circular arcs was 
attempted first. Its development is presented in Appendix 
B. It was anticipated that with the requirements of radial 
inlet flow and tangential exit outflow the double circular 
arc method could be used in the design of radial-outflow 
turning vanes. Once the method was developed and exercised 
it was clear that more geometrical constraints were required 
to produce a nozzle with an exit which was both nearly- 
tangential and convergent. 

The double circular arc program generated vanes which 
appeared to converge within the passage and then diverge 
towards the exit of the nozzle. Therefore, the flow would 
choke within the nozzle and the outlet flow conditions could 


not be predicted very well. 
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Modifications were made to the double circular arc 
design to achieve convergency at the exit. The first step 
was to remove the pressure line and use the camber line as 
the pressure surface. The assumption was made that the 
effects associated with inlet flow incidence resulting from 
the removal of the pressure line would be negligible. The 
results indicated that the passage was indeed converging 
toward the exit. The second step in an attempt to assure 
tangency and convergency at the exit was to allow for a 
thickness and a straight section at the leading and trailing 
edges. The straight section was such that a specified 
throat dimension was achieved at the trailing edge of the 
adjacent previous vane. The resulting program, which is 
presented in Appendix C, designs a wedge shaped arc blade 
such that a passage which converges at the exit is realized 
with a specified throat dimension. Using the method a high 
speed CDTD blading design was obtained. Figure 7 shows the 


resulting blade design. 
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IV. ELOW ANALYSIS 


A. METHOD 

A FORTRAN computer program for solving numerically the 
two-dimensional Euler equations using the Godunov method was 
developed at the Turbopropulsion Laboratory by Dr. 5S. 
Eidelman [Ref. 12]. With Dr. Eidelman’s help the prograa 
was to be applied to obtain a prediction of the flowfield 
through the vane passage. An abbreviated description will 
be given here of the computational method. 


The unsteady two-dimensional Euler equations are written 


as 
da + dOb+ B8c=O (7) 
Oc ox oy 
where 
e PU Pu 
a =/|/uUs, b= |p+ CUZ! , c =| PUV 
V UV p+ OP Vv2 
2 Ce+pdU Ce+pdV 
and 
2 2 
e= @e€ € + (U0 + V'D/2 > = energy of a unit of volume 
€ = p = internal energy 
CT “10 


€ = density 


U and V are the velocity components corresponding to the x 


and y coordinates 


U0 
il 


pressure 


~*~ 
li 


specific heat ratio. 
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Equation (7) is nondimensionalized using arbitrary 
reference quantities, and the Godunov method (Ref. 12] is 
applied to compute at each time interval the fluxes at the 
boundaries of the cells of a grid established for the 
problem under consideration. The steady state solution is 
found when no further change occurs between successive time 
steps. 

The computational domain is shown for a typical cascade 
problem in Figure 8 and the grid and cell notation are 
defined in Figure 9. Referring to Figure 8, for the 
tangency condition at the wall of the passage a solid wall 
boundary condition is imposed over the segments 3-5 and 4- 
2. Boundary conditions for segments 1-3, 5-7, 2-4 and 6-8 
are for periodicity of the flow with respect to the y axis. 

The parameters held constant at the inlet are the x 
components of velocity, U, the flow angle, 6, the gas 
stagnation enthalpy, H, and the entropy, S. The velocity, 
density and pressure can be calculated from these quantities 
easily as Vin» Pin, and pine At the outlet boundary only 
pressure is defined and other flow parameters are calculated 
in the solution. In the solution procedure, Uogyt, Vout, and 

wane are set equal to the values of U, V and eC one point 
ahead of the outlet boundary. 

A subroutine (program) constructs a non-orthogonal grid 
covering the computational domain as shown in Figure 9. The 


flow parameters related to the center of the cell have a 
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fractional index: j+%, k+%. The cell boundaries have one 
fractional and one integer index: j3, k+% or m+%,n. The 
parameters at the time t will have subscript indices and at 
time t+&t the parameters will have superscript indices. 

The Euler equation is approximated to first order by the 
equation in table 2. Solution to the Riemann problem is 
initiated assuming that the left state is that located at 
3-%,k+% and the right state is that at j+%,k+%. 

The solution to the Riemann problem will give 


(R,P,Un,Vt,? 


3,K+ % 

Un = the velocity component normal to 
the edge (3,k+1) -(€j3,k) 

Vt = the velocity component tangent to 


the edge (3,k+1) -(€j3,k) 
The capital letters R (rho), U, V, P are the parameters 
calculated at the edges of the zones. 
The component of the velocity is then transformed back 
to the cartesian coordinate system and a final value of the 
parameter on the cell edge is received 


(R»P,Un,B) yy 


B. RESULTS 
The program CEGE) was run on the IBM 370-3033 at the 
Naval Postgraduate School. Figure 10 illustrates the grid 


covering the passage for the wedge-arc blading shown in 
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Figure 7. An enlargement of the exit area is shown in 
Figure ll. 

In order to apply the code to the radial geometry, some 
modifications were necessary. A transformation of 
coordinates was made to conveniently accommodate the radial 
flow inlet boundary condition. Since the program remained 
basically two-dimensional, the outlet flow was treated as 
being periodic over the blade space at the exit which then 
did not change with radius. Thus the solution found will be 
valid for the flow within the rim and approximately true 
just outside the rim. For this reason the computational 
domain in Figure 11 was not extended radially outwards 
beyond one blade space. 

Results could not be obtained for the flow through the 
passage within the time available. Some difficulities were 
encountered by the program in the region of the throat. 
Some modifications of the grid must be made before a 


converged, steady flow solution can be obtained. 
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V. SUMMARY AND DISCUSSION 


Two exploratory projects at the Turbopropulsion 
Laboratory, 1). to test the concept of detonative combustion 
as an alternative to mechanical compression in gas turbines 
and 2). to design a teat rig for radial diffuaerasa uaed in 
centrifugal compressors, were examined. The projects were 
Similar in nature since they involved the generation of 
nearly-tangential radial outflow and each required a design of 
turning vanes for a test apparatus. The overriding limitation 
for either project was that any apparatus developed would have 
to be easily and inexpensively manufactured in-house. 

In the Rotary Detonative Wave Engine, “RDWE", project, 
turning vanes were required such that upon detonation the 
maximum possible torque impulse would be obtained. The best 
design possible to yield maximum power hence maximum torque 
was to produce turning vanes that exit nearly-tangential. 
Using Euler’s Turbine equation a relationship was developed 
relating exit angle to specific power and efficiency. 

The results for the unsteady expansion process showed 
that the value of the relative flow angle (82) established the 
limiting conversion efficiency and specific power. If the 
relative velocity increases beyond the limiting point the 
specific power generated will be absorbed rather than 


generated and efficiency will decrease. High efficiency 
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occurs in large part only when the specific power is low or at 
the limiting point. In general the conversion efficiency will 
vary during the expansion process, but will always be poor. 

This result implies that the concept of the RDWE can not 
be efficient as a stand alone rotor. Other methods would have 
to be employed to increase efficiency, such as the utilization 
of additional components, particularly those which would allow 
nearly constant flow by sequencing the detonation wave with 
respect to partial ports. 

The second project involved the development of a test rig 
for a high speed Centrifugal Diffuser Test Device, "CDTD". 

The project would closely parallel an in-house concept of a 
low speed CDTD. The low speed version simulates the flow of 
air delivered to a diffuser as a steady process. The high 
speed version simulates in a steady process, the transonic 
Mach numbers characteristic of advanced centrifugal 
compressors. 

For the high speed CDTD a passage would be required which 
would cause the flow to converge and choke at the exit of the 
passage. Because the concept of the RDWE was questionable the 
main interest tuned to designing vane shapes and the apparatus 
for a high speed CDTD. 

Two methods and their respective computer programs were 
developed. The first method involved the development of 
double circular arc vane shapes. Given the outer radius, the 


inner radius and the number of blades the radii of the suction 
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and pressure sides were determined. The design emphasis was 
on acheiving a tangent suction surface at the exit and a 
radial camber line at the inner radius. The concurrent 
condition of a tangent camber line at the exit gave the 
solution of a blade of zero thickness. 

The geometrical description resulting in a tangent camber 
line at the exit, resulted in the sine of the camber angle g@ 
being equal to the cosine of the subtended blade angle ¥ . 
Following manipulations to relate the value of #® to just the 
radius ratio (the only specified parameter), an equation of 
fourth order in sin F resulted. Using Newton’s method of 
interation the first (smaller) root was found. Using this 
value of ¥ a blade shape was generated. It was clear that 
the suction side was tangent but the efflux angle « was not 
zero. An efflux angle of zero would indicate that the camber 
line was also tangent to rim. A graph of the fourth order 
function indicated that two roots were possible for each 
radius ratio. Only the first root was solved for. It is 
probable that the second root would yield the tangency 
condition of both camber line and suction side. The selected 
root was shown however to give sin % = cos joe 

It was realized that convergence at the exit could not be 
guaranteed with the double circular arc blade design. 
Modifications were applied to the design to Paunaehe passage 
to converge. Improved convergence resulted when the original 


camber line and the suction line were used to describe the 
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blade. The simple double circular arc design required curved 
or thin trailing edges. With the addition of a finite 
thickness at the leading and trailing edge, a straight line 
segment, and a wedge angle all of which would better direct 
the flow in the tangent direction, practical blade shapes with 
convergent passage were acheived. 

Geoaetrical equations were developed to yield the proper 
blade design. The exit throat dimension was fixed in the 
design based on the limitations of the available air supply 
and based on dimensions scaled from a high speed CDTD design. 
It was anticipated that the flow would choke at the specified 
throat dimension between adjacent suction and pressure sides 
and the straight line segment would intersect the exit circle 
which was tangent to the suction side. This did not occur. 
The length of the straight line segment varied between 1 and 5 
percent of the value needed for this condition. It was not 
required that the exit circle be tangent to the pressure side. 

From the geometry and the specified parameters a blade 
was developed giving convergency at the exit. The camber 
angle % and the subtended blade angle ve were specified from 
experience obtained in the development of the double circular 
arc design. The chosen angles ¢ and ¥ gave an efflux angle ¢ 
between 8 and 93 degrees. To achieve the perfect condition of 
tangential flow angle « would be pacteead to be equal to 
zero. Based on current technology and manufacturing abilities 


an angle in the area of 8 or 3 degrees would give as close to 
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tangential flow as is possible. It is noted that the angle 32 
is related uniquely to the angle -. Once the radius ratio is 
given. 

The coordinates of the designed passage were then input 
into the FORTRAN program (EGE) to predict analytically the 
flowfield through the vane passage, by solving the Euler 
equations. The program (EGE) was designed to run on the IBM 
370-3033. The program successfully constructed a non- 
orthogonal grid covering the domain, however results could not 
be obtained for the flowfield through the passage in the 
initial attempts. It 
is thought that the difficulty lies in the curvature and 
density of the grid in the region of the exit throat. Since 
the EGE program has not previously been applied to a radial 
outflow problem with severe curvatures the task of analyzing 
the present flow with the EGE code is now thought to require 
the concentrated effort of the code’s author. 

In contrast, the CDTD vane design can be manufactured 
without difficulty and the flow from the turning vanes can be 
measured experimentally. It is suggested that such an 
experiment would provide a severe test of the ability of the 
EGE code to simulate transonic radial outflow containing 


embedded waves. 
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VI. CONCLUSIONS AND RECOMMENDATIONS 


A. CONCLUSIONS 

The following conclusions are offered: 

The Rotary Detonation Wave Energy concept in its original 
configuration is potentially inefficient. An increase in 
efficiency would require the addition of a component which 
could recover the unsteady kinetic energy of the flow leaving 
the rotor. 

A program for designing radial outflow vanes using double 
circular arc blades was demonstrated successfully. While the 
method may be useful in other contexts, it was found to be 
unable to produce converging passage shapes with a near- 
tangential outlet flow direction. 

A program for designing radial outflow vanes using wedge- 
arc blade shapes was used successfully to produce converging 
passage shapes with prescribed throat widths and with a near- 
tangential outlet flow direction. 

Both blade shapes can be produced relatively inexpensively 
from standard metal stock. 

A preliminary design of a test model of a high speed 


Centrifugal Diffuser Test Device was produced. 
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B. RECOMMENDATIONS 

The following recommendations and areas of further 
investigations are suggested: 

An arrangement of partial porting and following power 
turbine for the Rotary Detonation Wave Engine should be 
examined. 

The program for double circular arc blades should be 
modiied to seek both roots of the condition for camber line 
tangency, and to offer the selection of any value of ¥ between 
them. This will result in the ability to examine any value of 
Ww within the range of any practical interest. 

In the case of wedge-arc blades, the possibility of 
solving for the particular case in which the circle at the 
passage outlet is also tangent to the pressure side should be 
examined. A procedure should also be introduced to calculate 
the progression of the passage area taking into account the 
increasing radius. 

The proposed blading geometry should be constructed and 
the outlet flowfield examined experimently. 

Work should be continued to obtain a computational 
prediction of the flow through the passage. The model test 


could then serve as a verification of the code. 
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TABLE 1. Required Design Dimensions for Convergent 
Throat in High Speed CDTD. 





NUMBER OF BLADES THROAT DIMENSION SPACING 
Zz a Cin.) Ss Cin.) 
6 eS 4.1888 
12 ~-2613 2.0944 
18 -1742 1.3963 
36 -08708 -6981 
45 ~06967 25585 
60 05225 ~ 4189 
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TABLE 2. First Order Solution of 
Euler Equations 


a J+, k+ 5 ao. Ue sas At 
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DRT pt kh 


(B. ee 


jes j+y,k Sree k AX 


jtk,k 7 


AY C 


Beak kts ee ee ec) | 


C5455, k+1 Saye e | - Pith k+l : ee ce) v 
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Figure 4. CDTD Flow Generator - Section AA 





Figure@ #5. 


Qa 


Two-Dimensional Turbine Cascade. 
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Figure 6. Axial Turbine Blade Geometrical Design. 
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F= 9.1.72521902 ¢ = 80.82738092 


Condition of Tangency is 
sin $ = cos 

sin 6 = 0.159407755 cos z 0.159409430 
Chord # 0.983540343 chord test = 0.983540073 
alpha 2 18.34514098 

radial solidity = 1.909859317 

thickness to chord ratio ## 0.194933797 
radius of camber circle = 0.758551242 
radius of suction circle # 0.575175292 
radius of pressure circle # 1.462802226 


Figure 7. Final Blading Design. 





Figure 8. Computational Grid and Boundary of the 
Computational Domain. 
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a) 








b) Single Cell 


Figure 9. Computational Grid Notation. 
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Enlargement of the Grid in the 
Throat Region. 


Figure 11. 


APPENDIX A 
TORQUE PRODUCED BY DISTRIBUTED GAS REACTION 


A.1 METHOD 

The RDWE concept involves an unsteady process. The 
potential performance of the RDWE can be understood, 
approximately, by examining the range of conditions through 
which the unsteady process passes. Thus the steady flow 
performance of the proposed arrangement needs to be 
examined. We first attempt to predict analytically the 
power or torque impulse that can be generated from the 
engine under steady flow outlet conditions. 

The passages in the RDWE exit nearly tangentially 
(Figure 1). Since this condition is difficult to achieve we 
examine the effect of the exit angle on the amount of power 
that can be generated and the efficiency that can be 
produced. 

The basic equation for relating the fluid flow 
conditions to the energy transferred to the shaft is Euler’s 
Turbine equation. The torque of the fluid acting on the 
rotor is given 


T=m Eyev 97 - Fro V 02) AC1) 


If the rotor has a rotational speed w’, the shaft power 
transferred or rate of work transferred from the fluid to 


the rotor iss given by, 
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P =w’ T AC2) 


The work per unit mass flow is given by 


P/m =w’ (ry V0q1 - ro VOD) AC3) 


The subscript 1 denotes entrance station, 2 denotes exit. 
Defining the rim speed as U = rw’, equation AC3) becomes 


P / m = Wy = Uz VOq -U> Vo2 AC4) 


A general velocity diagram for a rotor is shown in Figure 
Al. 

The tangential component of velocity V9 is taken 
positive in the direction of the rim velocity U. Figure Al 
shows all positive velocity components in conventional 
cylindrical coordinates. 

The gas enters the wheel near the center, effectively 
with zero moment of momentum Cabout the axial direction). 
Therefore V61 = 0, and from equation A(4), 


Wx = P= -UD Ve ACS) 
m 


For Wy to be positive, it is seen that V60 is negative. The 
exit velocity diagram representing the flow at the rim of 
the wheel is given in Figure A2. 

From the geometry shown in Figure Al, in general 


sin 8 = V6 - JU AC6) 
iW 








Figure Al. Velocity Triangle of Conventional 
Cylindrical Flow. 





Figure A2. Velocity Diagram at Exit of Rim. 
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therefore 


Ve = W sin 8 + U A(C7) 


Substituting equations AC(7) into ACS) we obtain, 


Wy = = - U2 € W sin 8 + U Jj 


== 
Ta 


= ico CW sin § + 1 J 


U2 


i CW sin (-f - 1 J AC8) 


U2 


Defining the power conversion efficiency "c as the 
fraction of the available power actually converted to shaft 
work, thus 


No = P/m ACS) 


P/m + eee 


and defining 82’= -f2, the relative flow angle at the rim 


which is positive in the direction of rotation, then 


P = U2 { Wo sin 82’ - 11]. AC10) 
m U2 
From the geometry shown in Figure Al, 


Vo7 = Vieoe™ ¢ pe tv ope as 


so that, 


KVo+= (Up Vool= sUo- ae) AC11) 


Using equations ACS) and A(10) 
2 ; 
- U2 V82 = -Wy = U- C_Wo sin 82’ - 1] AC12) 
U2 


Substituting equation A(12) into A(11) gives 
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or 


zt 


By = % U5" [-Zeweesin io" + 1+ Ws e J AC13) 


U2 U2 


Substituting equation A(10) and AC13) into ACS) 


fc = P/m 
a 2 
P/n + VO 
~ U 2 CW. /U_) sin 8B.’ - 1) 
> ieee as 2 
U~ [(Wo/U2)> sin Bo’ - 11+%UD° [-2WsinBo’+1+ Wo “] 
U2 U2 


= (Wo/U>) sin 2’ - 1 
eC aon) 


Le - 2 € (Wo/U2 sin R2’ - 1 3 AC14) 


C 1 - W2/U2) 


Equations A(10) and A(C14) are valid for a discreet jet, 
for distributed gas at the rim. For the distributed gas, 


is possible to calculate the referred power as a function 


of wheel radius, giving specific power, Ps. Using equation 


AC10), 


Pa = P/mU2 = [(Wo/U>) sin B2’- 1] A(15) 


Therefore equation A(14) becomes, 


ic = econ ce 
(Wo/U2)2 - 1 


For sin 82’= 1, implying that 8>2’=90 degrees, 


c= 2 
(Wo/U2) +1 
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Thus the conversion efficiency is unity when Wo/U> = 1 
but this is a condition which, from equation A(15) gives 


zero power output. 


A.2 RESULTS 

The value of the specific power and conversion 
efficiency ‘(given by equation A(15) and A(16)) as a function 
of the ratio of relative flow to wheel velocity, for various 
values of the relative exit angle, are shown in Figure A3 
and Figure A¢ respectively. 

Figure A3S shows that the specific power increases 
linearly as the exit velocity increases, and decreases as 
the relative exit angle becomes less tangential. The 
conversion efficiency is seen to decrease as the power 
increases and to decrease as the relative flow angle 
decreases (less tangential). 

The results shown in Figure A3 and Figure A4 can be 
used to gauge the performance of the unsteady RDWE. In the 
unsteady expansion process the value of the relative flow 
angle (82) will be constant, being determined by the passage 
exit geometry. The wheel speed will be approximately 
constant and only the value of the relative velocity will 
change. Thus the unsteady process at different times will 
be described by points on a single line shown on Figure A3 


and A4, depending on the particular value of 82. 
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The difficulty, clearly, is that the exit velocity must 
not be allowed to decrease below a specific factor times the 
wheel speed or power will be absorbed rather than generated. 
Also the conversion efficiency will vary during the 
expansion, but will always be poor. 

This implies that the concept of the RDWE as shown in 
Figure 1 can not be efficient. The device would require 
some type of diffuser that would convert the kinetic energy 
leaving the rotor into static pressure. Unfortunately, a 
variation in the relative velocity magnitude implies 
variation in both magnitude and direction of the exit 
absolute velocity. A suitable diffuser for such a flow does 
not exist. Consequently, other geometries, particularly 
those which would allow constant flow to the diffuser by 
sequencing the detonation wave with (partial) ports at the 


rotor exit, should examined. 
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TABLE Al 


Specific Power and Power Conversion Efficiency 
Vs Velocity Ratio Data 


60 DEG 70 DEG 
LoS Pa Ne Wo/U>___—~Pa Ne 

0 4 2 0 -1 2 

55 -.56699 1.51200 5  -.53015 1.4137 
1.0  -.13397 © 1.0 -.06031 © 
1.5 .29904 .47846 1.5 .40953 65526 
2.0 .73205 .48803 2.0 87938 58626 
2.5 1.16506 .44383 2.5 1.34923 .51399 
3.5 2.30109 .36108 3.5 2.28893 .40692 
4.5 2.89712 .30099 4.5 3.22862 33544 
5.5 3.76314 .25731 5.5 4.16831 .28501 
6.5 4.62917 .22444 6.5 5.10800 24766 
7.5 5.49519 .19892 7.5 6.04770 21892 

80 DEG | 30 DEG 

/U> Pe Ne Wo/U> __—~Pg Ne 

O -1 2 Oo -1 2 

5 -.50760 1.35359 sas 1.33333 
1.0 -.01519 © 1.0 O 1 
1.5 .47721 .76354 1.5 5 . 80000 
2.0 .96961 .64641 2.0 1.0 66668 
2.5 1.46202 .55695 2.5 1.5 .57142 
3.5 2.44683 43499 3.5 2.5 44445 
4.5 3.43164 35653 4.5 3.5 36364 
5.5 4.41644 30198 5.5 4.5 .30769 
6.5 5.40125 26188 6.5 5.5 . 26667 
715 @enaeeae 23117 7.5 6.5 23529 
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APPENDIX B 
DOUBLE CIRCULAR ARC BLADING DESIGN 


B.1 METHOD 

A computer procedure was developed to design radial 
out-flow, near tangential turning vanes, specifying only the 
outer radius, the inner radius, and the number of blades. 
The blade spacing, radial solidity, camber angle, suction 
and pressure surface radii, chord length, thickness to chord 
ratio, and efflux angle are determined. The suction surface 
is designed to be tangential at the rim. 

A Hewlett Packard 9830A computer was used. The blade 
shapes were generated using a Hewlett Packard 9862A plotter. 

The geometrical approach to the design of the vanes is 
illustrated in Figure Bl and much of the notation is defined 
in this figure. The following assumptions were made as 
initial constraints: 

1. radius at exit Ro, is known 

2. radius at entrance Ry, is known 

3. number of blades Z, is known 

4. The suction side is tangent to the outer 

radius (point C). 
3. The camber line c’, is tangential at the inlet 
radius (Cpoint A). 


6. The blade subtends an angle VW at the center. 
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Initially, a procedure was developed which required 
specifying the camber angle %, Ro, Ri, in addition to the 
number of blades. Since the choice of % which would give 
acceptable values of efflux angle and blade thickness was 
not obvious, an equation for the condition that the camber 
line was tangent at the rim (which would give a blade of 
zero thickness, since the suction surface was also required 
to be tangent) was written. The equation was seen to have 
two possible roots. The selection of the smaller 
Cincorrect) root gave positive values of the efflux angle 
and reasonable values of Mitade thaenesee The program which 
was developed provides the option of selecting @, or of 
selecting the blade which results from finding the incorrect 
smaller root to the tangent camber line condition. 
Described here are the equations for the generation of the 
blade shapes. 

Defining peripheral blade spacing s as, 

& = 21 Ro / Z 
and the angle subtended by the blade space as 
Ys = 3260/2 
the radial chord as 
C= (Ro - Ri? / Ro 
then the radial solidity is given by 


Sr = (Ro - RjDO / S&S = 7 Y 
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Using these expressions and Figure Bi the following 
equations can be written: 


GE = Re - Re cos ( £/2) BC(1) 


Re (1 - cos ( @/2) 
HE =) Gaiw2 + Re-t 1 = cos: 4/21 BCZ) 
Where Y is the thickness to chord ratio. Also 
HE = Re € 1 - cos (#8/2)] B(3) 
To solve for the thickness and YT we subtract the 


equations B(1) and BC2) 


HE - GE = t/2 = ¥ c/2 B(4) 

Using Figure Bl, 
Pe hel CoS © S/2)) = + c/2 BCD 
= Rp (1 - cos (%p/2)] BC6) 


Subtracting equation BCS) from B(1), 


GE bee t7205= 7% c/2 BC7) 


The chord c, is given by, 


c = 2 Rp sin (%p/2) B¢8) 
= 2 Rs sin ( 4s/2) B93) 
= 2 Re sin (8 /2) B(10) 


All equations were made dimensionless by dividing by Ro, 


giving for example 


c/Ro = c" B¢«11) 

c” = 2 Rp sin (#p/2) BC12) 
Ro 

= 2 Rea sin (95/2) BC13) 
Ro 
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= 20R5 ein Go B(14) 
Ro 


For the geometry we seek, given Ro, and the camber 
angle %@ and efflux angle a, the angle ACO (or ) is given by 
Nn =90 - €¢ X£ + S/2) B15) 
From the triangle ACO 


y+ 90 - (a + 9/2) + 180 - g/2 = 180 


resulting in, 


90 +F = G+ B(16) 


This equation fully defines «. The values of % and &% 
are needed to complete the description of the camber line. 
To solve for the radius of each blade arc the following 
relations were used. 


First with 


As = Re - Re cos £/2 = HE B(17) 
Ro Ro 
and 
Ap = Rp - Rp cos (%p/2) = FE B(18) 
Ro Ro 


Using equation BC(3), BC(6&) and BC7) 


Re € 1 - cos (8/2)] - Rp € 1l-cos ¢(%p/2)] = ¥ ¢c/2 B¢(19) 


Also, 


Ag - tT c/2 - Rp €1 - cos (%p/2)] =  c/2 B(20) 
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which gives 
hPa Aa YC B21) 
Using this relationship between Ag and Ap we can solve for 
Ra and Rp. From equation B(2), and BC(3), 
Ag = Rs [1 - cos (%s/2) B<2Z2Z) 
or 
cos (%s/2) = 1 - As/Rgs 
Squaring and adding together equations B(22), and B(10) 


we obtain 


cos (gs/2) + ae (85/2) = (1 - ee + (e/ PRs) 
or 
1 = 1 - 2¢€As/Rg) + (Ag/Rg)2 + C2/4Rse 


which becomes 





Re = “BMGgee oc * 1 B(23) 
BA. 
or 
2 
Rei = WAm + ¢c ~ 1 B(24) 
Ro Ro 8As 


Using equations B(S), BC(6), and B(16) a similar manipulation 
will yield, 


Rp = %* Ap + <7 i B(25) 


The angles %s and %p are found using equations B(13), 
and B(14) since 


Sin (9s/2) = c/2Rs 
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or 


Gs = 2 arc sin (¢( c/RO) B(26) 
2Ra/Ro 


and similarly 


Sin (Sp/2) = c/2Rp 


or 
Gp = 2 arc sin (_c/Ro) BC27) 
2Rp/Ro 
The desired condition is that the suction surface be 
tangent to the outer rim. From Figure Bl this is given by 


. the relation 


Gs/2 = X~ + G/2 BC28) 
Y is not an independent parameter since it can be 


calculated using equation BC(25), thus 


y¥ = C (Rea C€ 1 - cos (%s/2)] - Re € 1 - cos (%/2)] ) BC29) 
(c/2) 


To insure near tangential flow at the exit we first 
seek the condition of a tangent camber line at the exit. 
This implies since we also seek a tangent suction surface, 
that the blade would have zero thickness. Using the 
notation in Figure B2 the following equations apply: 

Sin (8/2) = L/c BC30) 
and 

sin ¥ = L/ Ro B(31) 
Substituting equation BC(30) into BC(31) gives, 


c = Re sin ¥ B(32) 
Sin (8/2) 
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From the geometry 


cos (8/2) = y/e Bicas) 
and 
cos ~ =_Ri + y B(34) 
Ro 
or 
Ro cos ® = RG + y B(35) 


Substituting equation B(33) into BC35) gives, 

Ro cos # = Ry + ¢ cos (8/2) B(36) 
or 

Ro cos - - Ri = ¢ cos (6/2) B(37) 
Equation B(32) can be rewritten such that, 


Ro sin ¥ = ¢ sin (8/2) B(38) 


and dividing equation B(38) by B(39), 


Ro sin ¢ = tan (6/2) BC39) 
Ro cos ~ - Ri 


This equation establishes the relationship between ~ 


and #. Using Figure B3 the relationship of ec, ¥, @ is 


established, 
cos (8/2) = _¢c/2 B(40) 
R tan vw 
or 
c = 2 Ro tan ie cos (6/2) B(41) 


This is the condition for the camber line to be tangent at 


the exit. 


Using equations B(32) and B(41), 
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sin ¥ = 2tan # cos (9/2) B(42) 
sin( 6/2) 


giving, 


cos ¥ 2 Bin (8/2) cos (9/2) 


2C % sin ¢€ $/2 + $/2) + Sin (8/2 - $/2) J 


sin B(43) 
Therefore for the camber line tangency condition 


cos ¥ = sin @ B(44) 


Using equation BC(39) and B(44) we obtain the following 


2 sin ($/2) cos (8/2) = tan (4/2) 
2 cos 2 (/2) 
and 
2 sin (%/2) cos (#/2) = Ro sin ¥ 
2 cos 2 ($/2) Ro cos ¢ - Ri 


Using the general trignometrical relationship, 


Sin 2x = 2 sin x cos x 


we can manipulate the preceeding equation to obtain 


Sin % C€ Ro cos ¥ - Rj; J] = 2 cos 2 (8/2) B(45) 


Using equation B(43) this becomes 
cos ¥ [€ Ro cos ¥ - Ry 1] = 2 cos@ (9/2) B(46) 
Using the relationship 
cos? x = sin? x + cos 2x 


2 cos2 (/2) 


2 (Csin2(%/2) + cos 9] 


2 €_1 + cos 9% Jj 
a 
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= 1+ (1- sind a) % 


Using equation B(43) this becomes, 


2 cos2 (¢/2) 1+ (i =- cos2s) * 


1+ sin ¥ B(47) 


This gives 


Ro cos* %¥ - Ry cos ¥ = Ro sin # + Ro sin2y B( 48) 


Writing s=sin WY equation (48) becomes, 


Ro(1 - s2) - Ry ¢ 1-82)” = Ro 3 + Ro ac B(49) 


or 


Ro (€ (1-82) - s2 - 8] = 452) 


Ry 


(Ro)2 € 1 - 28 - 382 + 4683 + 484 ] = 1-52 
(Ry) 


Thus equation B(48) gives the condition 
Fis) = (Rj)2 (1-82) - 1 + 38(2 + 3s - 482 - 483 ) = 0 B(SO) 
Ro 

The condition of tangency of the camber line is given 
by the solution of equation B(S0) for s between O and 1, 
since s=sin ~ . Using Newton’s method of interation the 
root can be found for a given radius ratio. Figure B4d shows 
FCs) in equation B(S0O) for various radius ratios. Figure B4 
clearly identifies two possible roots. The larger root is 
the correct root, giving zero blade thickness. The smaller 
root if incornacty but gives positive values of G and finite 


blade thickness. Emphasis was placed on finding the first 
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root. The value found will give the selected value of sinf. 
Once * is found, using equation B(39) g can be calculated 
and c is obtained using equation B(32). 

Using equation B(28) and B(9) gs and Rs can be 
calculated and the thickness is then given using equation 
B(29). 

A test is made to assure that the length Ap is 
positive. If Ap is negative this implies a convex rather 
than concave pressure side. This would also imply that the 
inner radius is too small and should be increased to produce 
concave curvature on the pressure side. Ap is found, from 
equations B(19) and B(16). 

The alternate procedure did not require the camber line 
to be tangent at the exit. In this case the value of ¥ is 
not solved for but is left as an input parameter. The 
values of # and c are then found directly using equation 
Bsa). 


The centers of curvature are given by, 


Xds =_Re - (_Re cos (8/2) -_Re cos (%s/2) ) cos (8/2) 


Ro Ro Ro 

Yds =_Ri + (_Re cos (8/2) - Re cos (45/2) > sin (8/2) 
Ro Ro Ro 

XqdP = Re + (_Rp cos (%p/2) - Re cos (8/2) ) cos (6/2) 
Ro Ro Ro 

YdP = Ri - ¢ Rp cos (%p/2) - Re cos (#/2) ) sin (8/2) 


Ro Ro Ro 


with 
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Xq = Re 

= oRe/ Re 
mas Ri 

= Ry/Ro 


To draw successive blades on a rotor the program must 


rotate the axes about the rotor center after each Dlade is 











drawn. The following equations were used (see Figure BS). 
X = x + Y sin 9 
cos @6 cos @6 
Y + x Bin 9 = y 
cos 06 cos 86 
therefore, 
X = x cos ®@ + y sin 6 
Y = y cos 869 - x sin 8 


Given the number of blades, the following operations 
(for all three arcs) are performed, the coordinates are 


calculated as 


x = Xd»ds»dP - R cos 6 


Y = Ydrdsreds + FR sin 6 
The arcs are then drawn from 6 = ( (8/2) - (%s,p/2) ) to 


(CG/2) + (88,p/2)) 


using 
X = x cos (N) 360 + y sin (N) 360 
Z Z 
and 
Y = y cos (N) 360 - x sin (N) 360 


a2 ie 


where N is the number of blades. Figure B6 and B7 
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illustrates the computer generated shapes and passage for a 


60 blade design. 


B.2 COMPUTER PROGRAM “ACIRC” 


B.2.1 Description 

A computer program, “ACIRC” was written to 
construct double circular arc blading using the above 
procedure. The program will calculate and draw blading for 
@ initially radial-outflow using Newtons Method to solve for 
@ suitable value of ¥, or for an arbitrary input value of ¥. 
In either case the method will yield a duct 1sawlemeame which 
is tangent to the outer rim. 

The program is user friendly and is annotated for 
ease of understanding. The program is divided into six 
major parts. 

The first part is the input data section. The 
second part calculates the periodicity, radial solidity, 
spacing, and asks the user if the value of ~ is to be 
calculated or to be specified. 

The third part calculates the first root for 
tangency conditions at the exit. No input is required for 
this section. The fourth part calculates the parameters 
needed for the suction surface to be tangent to the outer 
rim. Usually, if the user specified a value of # there is 
the possiblity of the distance from the chord to the 


pressure circle being negative. If this occurs the result 
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is a convex arc. To remedy this, the program asks the user 
to increase the dimension on the inner radius by a 
percentage. The program will use this value to solve for 
the required inner radius to achieve a suction surface that 
is tangent to the outer rim. 

The fifth part prints the output. The sixth part 


draws the blade surfaces. 


B.2.2 Required Inputs 


The following inputs are requested: 
Outer radius 
Inner radius 
Total number of Blades 
Yes or No response for solving for the value of ¥ 
Percentage increase in the inner radius - may 
be required, when the value of * is specified 


Number of Blades requested to be drawn 


B.2.3 Listed Outputs 


The following outputs are listed: 
Outer radius 

Inner radius 

Number of Blades 

Angle e 

Angle @ 


Results of test for tangency 


el 


Chord length 

Angle ¢« 

Thickness to chord ratio 
Radius of camber line 
Radius of suction line 


Radius of pressure line 
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B.2.%4 Program Listing 


Listed here is the program listing,’ ACIRC” for 


design of double circular arc blades 


PEMSss*4422HC IFC 1 JUHE {asdeeeeex424% 
Beers ee ee PUL ie CAREUERIG & DRAW RRC ELADING++*+++<* 
REM+ex+4eeF0R AW THITIALLY OUTFLOW CASCRIE++s+++ 
REM*+*#+e22USITHG HENTOHS METHUD TO SULVE FUR THE+*#+++*+#+* 
RPEMSes20F PSl+sx+eeeee% 
DEG 

mbRee =Lals—-isl 
mise “SET FLUTTER 
STOP 
RPEMS<xS42THP UT DARTH ee ee ee te eee ee ee Se ESS SHE 
eee CHUERT IRE WUTER KHEIUS ; 

IHFUT RG 

HISP “ENTER THE IHHER RAUIUS’s 

IHFUT Ril 

R4=R1°R6 

DISP “EWTER THE HUMEER OF BLADES"; 

tae) 28 

= i 

BeMmeteeeeeeEHD IHNFUT DATRHsee% 
REMe*#4*CRHLCULAT IONS 


7 f- 


aon SEG a VEU) 7s 


eee SK He HHH ee EH eH HHH HH 


mn: 


Se.’ 


ut 


I= seh Zk 

pase LAL SS ees4 

S2=eePleRu- cu 

Be (Ree k iose 
Se=C(RU-R1d7Se2 
SW=SE°RB 

eS SUES] 
DISP “WAHT FSI CALCULATED Y=1,0=H"s 
WAIT 108 

ito Ot 


IF O1=1 THEH 348 

DISP “IHMFUT FSIS 
IHFUT 33 

O2Z=COSS3 
TO=ATHCOZ7S0R 61-0212) 
GOTO 330 

GOSUB 1358 
RPEMS*2+4¢HOM DETERMINE THE WALUES FOR FOIsFHIsAND C+#+%# 
SS=ATHCKS“SORCL-at2y) 


) 


PSone tolls s-OR- CBS sa—F 1) ) 
we= 5 1h5-3 
O6=COSTS 


Ci=RHs5SITIHSS“SIHCTO“2) 
RETIRE ERK E ERE FERS HES FEARS REESE EFS EEE EK KEHES 
RPEM+*4+**4£4TEST FOR TAHGENCY*+2+% 
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Ce=feTAHss+C0st ToS skin 

REMSS#e444C0HCLUIES CALCULATIUWH FUR PSIsPHI:s AH Cx+eeeee 
Ai=3 H+ s-Te 

eR 
M4“ CSSSTHC Tees 3) 

Ri-Ry 

DE=ATHRS“R AD 

PEPER H SHR SSSR SSE SHS SSS HSH SH eee eee eee 
REMS++sFO0R SUCTIUH SURFACE TANGENT TO ULITER FRIM+++# 
Pe=C2eno ie 

RS SE ee Sl ieee 

Rebterbie oe Ter 2 a ee eae 
ALig=RSec P-COos8 Tees oo4+lese Tle 
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TABLE Bi. 
for Double Circular Arc Blades 


*ACIRC” 


AO 


Al 


A2 


CO 


C1 


C2 


Ea 


2 


I3 


Ki 


K2 


K9 


NO 


N9 


Ol 


TEXT 


List of Computer and Text Symbols 


IDENTIFICATION 


efflux angle 

suction line to chord length 
pressure line to chord length 
c/Ro dimensionless parameter 
for chord length 

chord length 

chord length for tangent test 
Interval (8/2 - %s/2) 
Interval (8/2 + %s/2) 


Interval (98/2 + @s/2) - (8/2 
- $s/2) 


constant value 
Newtons method, y/y’ 


initial guess for Newton’s 
method 


counter 
Number of blades to draw 


Number of blades to draw for 
loop 


Number of points for interval 
of arc 


Equals 1 for solution of ¥,2 
for input ¥ 


TABLE Bl. 


“ACTIRC* 


O2 


O03 


04 


OS 


06 


Pil 


P4 


RO 


Ri 


R2 


R3 


R4 


RS 


R6 


SO 


Si 


S2 


S3 


SS 


TO 


List of Computer and Text Symbols for 


Double Circular Arc Blades (con’t) 


TEXa 


cos ¥ 
Re 
Ra 
sin 


cos @ 


Ri 


Re/Ro 
Ri/Ro 
Rs/Ro 
Rp/Ro 


Vee 


Or 


1) 


IDENTIFICATION 


cos ¥ 

output print of Re 

output print of Rs 

for tangent test 

for tangent test 

n 

constant, cRo/2Rp 

outer radius 

inner radius 

density, (Reo - Ry)/Ro 

camber radius, dimensionless 
inner radius, dimensionless 
suction radius,dimensionless 
pressure radius,dimensionless 
periodicity, 2MRo/Z 

radial solidity 

spacing 


blade curvature with respect 
to camber line 


periodicity, 360/2 


camber angle 


TABLE Bl. 


“ACIRC” 


Tl 


T2 


T3 


T4 


T9 


X1 


X2 


X3 


X4 


XO 


ved 


Y2 


Y3 


Y4 


Y9 


zO 


List of Computer and Text Symbols for 


Double Circular Arc Blades (con’t) 


TEXT 


gs 


Xds 


*xdP 


Xd 


Yds 


YdaP 


Yd 


IDENTIFICATION 


thickness to chord ratio 
suction angle 

pressure angle 

loop parameter to draw arcs 


tolerance on Newtons method 
of integration 


Suction center x coordinate 
pressure center x coordinate 
camber center x coordinate 


variable for x values in 
continuous loop 


variable for x values in 
rotational loop 


suction center y coordinate 
pressure center y coordinate 
camber center y coordinate 


variable for y values in 
continuous loop 


variables for y values in 
rotational loop 


number of blades 





X 


Figure B1. Double Circular Arc Blade 
Geometrical Description 
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Figure B2. 
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Figure B3. Continued Geometrical Description for Camber 
Line Tangency Condition 


82 





‘ 
mee- 7 ow 





i 





Se wea we 20) Ol eae) e 





en ee 


- —+ ——- 
pre be te 
=) 1G = =e 
Pe ee 





Se ee 





-o—- 4 - ee 


ap -~- ef wes wees wm mmm oe - - he 8 mee He - 8 = oe 





o— 8 epee eee 








- oe: 














a eo2e+- +2 = . - - -—- —~—e- - - -_—- 4 co ; ; = . 

hein ee Ce cane ERE. SE eee ‘ heir 

ent ae cree a8 - ns = Z 
ele ar a eg oes oO 4 

tea a. a ne es el cer eae . : 

a ens ; : : 

e2-- ne == - = e = - . é 









- ce & care 5 : rhagm _ 4 - at - os -¢ - - s 
— + ae eee - - - fF - -- - -_—— awa we - - 4- -2 2 ee - - = -=¢@--- °- 

a -- - $ ee eee - we hee - ~ 4 

a “- 2. (-- - 77 oe _— - tee -—— e 

— ee —_ 

- - ~ - = 6) ==. = oe a oe — - 

can om -« ere -- i . 4 

es es ee . 

-.- _- —_— —e - 6 - . 




















: . i ee - ¢ - : ) - - - e -_=— s “= e 
~_—- one ee ene Seo = oo === ee ee a ee ae Sa alae tip are A oe ee me -¢-3* Gh aa -° 1 = = 
« t0 4 seucé=« eee ew eee . SSS a. ee os a arse eh ae al ar oak ae are « eae tue owe 8 «= « - — - « iomeien - = how ee - e—<- 6 i eat, a) oe $ . --—- 2 -- -¢ 
ose 6 oe sb om ee tee ce ee = Rt ee mm hee ee eae ee oe eee «= 0 6 -- ae ae ya ~ 
ee ee -pere- = -- tec nee mis ome mS = One ‘ See ° —- aes e 2 . . — ae - . - - 
Dee ceo hat eec te oe Piel alot cere he ele tete wie = ea eet) ee = = — s J See 5 ee eee -- 2 t = = { Se . 
-- eee - pe -- ¢—— aye) velo too. ew ee ee - = 6 me = e -<2 = & - . --- - --- - --+e — - - 4 ---= a --- . 
- _ - }- a —_ -_—- - =: - Oe & =“. - e -- . - - - + & - --t - - —--- . - s 
.-- + - t- cee Ke = - @--- -- - e - - - - ' - . * =e: - 
aoe, ate ge = eee eee =: eeu ee ene ¢ -sm-- = -4 io <n re De oyie ! . 
a -- eee a —--e ee mee ee -- - -- — -*e = -e - - ra -- 1 - t ° - - + . - - - 6 
= mal een eas ’ Nae eon ce atte * = ‘ See 4=- a = ¢ e 6 e Biav ear se 8) mate Be ee Be 5 = ete terete a . - e 
a ee he i 
a ’ = - ¢t - --4 - ° ? - . - “8 . -_- + 
=e = 4—= . er es ett eee a pene ® eS -- . i ° 
o- 4 -- e- e 2-5 = - -- - .- -——_— - Siceyhe ela is ome ca melee : ® = = -_- = -- -~ 
ee - Ff ee ~— o- — --= e - o-—-- 4 -- - © . ee en 2S ® - - ses -—— - - - - 
LE 
ernme Se ee - cs 7 —< Sets :, a = t rte 4 = - Sn als ae e 
ees e ee a eT - pie eL -@se ee ee - * - - ' . - : -—=— s 
Sg Ce ee Cee eae 4 eee acwer/o ham Ae aos om--— + ft -- Sas ee -~ - . oo eee en e« — § --~+e----2 
- da-- eran mere —e- - wey oa Oe ‘ee > oe --- ee ---@- —-- 4 
ea a A 








« -¢-~- - em -e = ---e come = . 
e-- p= - eee . = i eee =e a a ee a 
mee pm em foe ner ee ee - - --2 = + - --- _—+ . 
: e- = en | ' oe C=) ona - -%- a eo . 
-- ~ a ° oe eee aie ° ip a --- ¢ 5 a a © 
- ne oat et eon 8 4 oe - { =m~ 2-2 —- - . — se ew eee eh me —- - 
-_— is - wmonm tea. m--- - - ms own - ee mw ee eh eee ew ele 

ores - ~ " ea --e7 = - ee -_eme er ewe ? 





e+e @ -- 











-=- OS OSS oO fp Se OG cme ne tee fF - e ee ere - 

-- —_—- - = te meme - = Fe ee Ke ee Oe = enone me we = -- eer =e + o-8 ee ? 
we ee me eh ee ee eee bw ese ee pe ee oor «we 4- « ? 
Bias eos ae ones ae ee ee hy et ees) ene eas wes, oe oe aie elle '= es 














o-oo He He ee aan —2 = © 2 owe ra cee eet my 4 a : 
ee Oe ee ee wom be 2 = = a ee eee QO - a- sn eee . ee - -—- + - oe mee § o- — eee ee - —-e- - - ° 
net —-= ee vara paee taney pea = = § ee ee -- —- Cetin = twit wd inn =~ = > woe hee ee Fe “- --- - 4 
ee -- a ee one - 
—-— -~ _-— -_ = -- mee a a ——.- ooo wee -- -_efe “ese 2s ee . -- ~ --2* t- - 2 
+ . ~—- ~~ — + -e Ww 4 + + = 
--¢- —-~ —-- - 4 aa ee cs 
Laka ce sib anes Dace |: at ee ee eee . 
ee te . ene teen ad) saci oe GS ee ee a eee ant aie ila one = 





- &-—«. ~ oo oO ee 


See ey el hee ee eee 





ee ee ee 





— Se - mee ee 
a 


Spee eeen eis whe) Pac sees eae 


-— cee Om -oe 





- mm - e-em -- ee we Hee oe 







Function of Y Resulting from the Condition 


of Tangency of the Camber Line (equation B(SO)) 


Figure B4. 
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Figure BS. Schematic for Rotation of Axes. 
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-® 


Ro =4 Ry = 3.2 No. Blades = 60 
Y= 9.1.72521902 £4 = 80.82738092 


Condition of Tangency is 
sin © = cos 

Bin 6 = 0.159407755 cos z 0.1359409430 
chord = 0.983540343 chord test = 0.983540073 
alpha = 18.34514098 

radial solidity = 1.909859317 

thickness to chord ratio = 0.194933797 
radius of camber circle = 0.7585951242 
radius of auction circle = 0.57535175292 
radius of pressure circle = 1.462802226 


Figure B6. Computer Generated Shapes 
for a 60 Blade Design. 
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APPENDIX C 
WEDGE-ARC BLADING DESIGN 


C.1 METHOD 

Because of the limitations seen in the outlet passage 
shapes which could be generated using purely double circular 
arc blades, symmetrical linear surface segments were added 
in an attempt to achieve a convergent nozzle at a shallow 
angle at the rim exit. 

From the design using only double circular arcs it was 
found that by using the original camber line and the suction 
line to describe the blade, improved convergency at the exit 
was obtained (see figure C1). Removing the original 
pressure line would cause the flow to be at a small 
incidence angle at the inlet, but this was thought to be 
acceptable. 

By allowing the flow to exit over a straight suction 
Side from the exit of the passage, better control of the 
outlet flow angle was assured. The wedge arc geometry is 
shown in figure C2. 

Referring to figure C2, the wedge arc blade shape was 
developed to have the following features, or constraints: 

1. The pressure side is a circular arc. 
2. The suction side between point H and I is a circular 


arc, between straight and symmetrical line segments. 
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3. A wedge angle (¢€) is specified at the leading 
and trailing edges. 

4. Equal leading edge and trailing edge thicknesse (5) 
are specified. 

S. The radius to the suction surface at the point of 
tangency to the straight segment, (point I) 
intersects the pressure surface of the adjacent blade 
at the trailing edge (point C’), or downstream of 
point C’, depending on the value of the parameter 


Craw 


Referring to figure C2 and C3 the following 
relationships are obtained: 
XF = X¥c ~ 8% cos @ C¢1) 


YF = Yc + 8’ sin @ C (2) 


where §’ is the leading edge and trailing edge thickness. 
The value of 6’ has been made dimensionless by dividing by 
Ro- All other quantities have been made dimensionless in 


the procedures given in appendix B. 


Xt = Xc ~ 68/2 cos @& CCSD 
Yt = Ye * 8/2 sin C4) 
Xe = ~§ C¢(S) 

Ye = Ri C(6) 

xL = ~8&/2 C(7) 

YL = Rj C8) 
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The wedge angle, € is specified. With € known, the 
length of the straight section (£) from the exit can be 
found. For a symmetric blade, the length (££) is the same 
at the inlet and the exit. Referring to figure C2, C3 and 


C4 

Q= Pea = €7) C9) 
where f = 0 ’ when the point I lies on the radius C’D’ 
(figure C4). The slope of the line C’F (Cy) is given by, 


y= arc tan (¢ yF - yo’) C¢c10) 
( XF - Xe") 





The coordinates of C’ are found by solving for the 
previous blade and rotating the coordinates. These values 
are placed in an array until needed. A straight line is 
drawn from the current blade to the previous blade as is 


defined as. 


C’ F saat ye yeti + (€xF - xo) ) * Cc11) 
Referring to figure C3, the slope of FG is given by, 
Azta- #@- €3) Gic1 2) 


A relationship of £7’, C’F and the line FG is given by, 
f’ = C’F cos C’FG 


= C’F cos € a - vad - € + pj Cc(13) 


From equation C(13) the coordinates of I and H are found, 


XJ = XE - L cos (a - e - €7 C(14) 
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YI = yr - @ sin fa- ¥ - €] C(15) 


xH=- €§ + £ sin 61 C(16) 


yH = Ri/Ro + & cos € COZ) 


The centers of the suction side are given by either, 


XD’ = XH + Re cos ¢ C(18) 
yD’ = yH + Rz sin ¢ _ 6619) 
Or, 
Xp’ = xT + Re sin C€ CG - em €] C620) 
yD’ = yr - Rez cos € a@- #% - €)7 C¢(21) 
From equations C(17), C¢(19) and C(11) R2 can be found 
since 
XI - XH = Rez (cos ¢ - sina ) 
yielding, 
Rz = MTs C(22) 
cos ¢ - sind 


To calculate the position of I for a given throat 
dimension (a’) at the exit, further manipulation is 
required. Referring to figure C4, where “a™ denotes the 
center of the exit throat, 

Xa = Xc’ + a’ Sin « C¢(23) 
and 
Ya = Yc’ - @’ cos ¢ C(24) 


From equations C(21), C(13) and C(15), 


Ro = ( xF - Q cos 4+ 5 + sin be 
(cos € - sin AD 
= xpo+ § + {sin € - cos4] C¢(25) 


[cos ¢(-sin qa) 
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Referring to figure C4 
C Re + a’J2 = € xp’e- xa 2 + FC ya - ypo’)2 C(26) 
Using equation C(20) and C(21) the above equation becomes, 
( Ro + a’)2 = (xz-Rz sin QA - Xal2 + Cyag-yr-Rz cos A)2 
referring to equations C(1i4), and C(15), 
[ Re + a’)2 = txp - fcos a + Re sin A -xal2 
+ € ya - YF * Asin a - Rz cos A j2 
which reduces to 
2Rza’ + a’2 = (xp - xQ)2 + (yg - ye2 + 2 C(27) 
+2 q cos A (Xq - xf) + 2 J Sin A (ya - YF? 
+2Rz cos A ‘(Ya - YF + 2Rz sin A (xF - Xa). 


Using equation C(25) the following results can be 


obtained 
a’2 = (xF- X_)2 +lya- yr] +Q2 + 2 cos A (xa-xF) 
+2 £ sin Alyg -yP? 
+ 2(xp_+ § + {tsin ¢€ - sos 41 ) sin A (xF-xa> 
[cos € - sin A J] 
+12 oth o2 6. + aan Coa cos 4 | | cos 4 (ya-yP? 
{cos € - sin 4 J 
+ gid + § + ian € - cos &) ) a’ 
{C cos €¢ - sinAJ 
[cos € - sin A Ja’2 = (€ cos € - sin QJ * C(28) 


(XP - Xa)2 + Cya - yRoe +K2 


+ Be eon ra xe oe 2) eindChecyE) 


Xf + §& + ean € - cosAjJ sinASC(xp-xa) 
xF + § + M€sin € - cosdA] cos A Cya-yF) 


+ + 


NN 


+ 2 xF + § + Nivea € - cosAJj a’ 
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and on combining terms, 


[cos ¢€ - sin i jaa = € cos € - sin A ] CeZzoD 
+ {2 (€cos ¢ - sin A J cos A (xa - xP) 
+2 (cos € - sinAJj sinAlye - yF) yX 
+ 2 xp Sin ACKF - Xq@g) + 2 § SinA (xp -xa) 
+ 2 {cos € - sinAl] sin & (xp - xa 
+ 2 xp cos Alya -yF) + 2 & cos ACyae - yP) 


+ 2 {cos €¢ - sinA]J] cosAClya - yP) 


* 2xp a’ + 2 § a’ + 2 8 [eos € - singe 
{cos €¢ - sin ] a’2 = € cos € - sind] 2 C¢(30) 
+ { 2 Ccos € -sin A] cos A (xe - xP) 


+ 2 €cos € - sinA] sind Cya- yP) 
+ 2 (cos € - sind] sin ACxF- xa) 
+ 2 {cos € - sinA] cos Alyse -yP) 
+ 2a’ (€ cos € - sind] } 
+ 2xp sinAC(xE - xq) + 2 § sin ACxF - Xa) 
+ 2 xr cos Alya - yFD + 2 § cos AlCya - yp) 
+ 2 ME VS + 2 baa 
The above equation is quadratic in the unknown g, and a 
solution can be written explicitly. The value of £ can be 


found using equation C(9). The result is as follows, 


Pa 


B= { 2 {€cos € - sinA] cos A(xg - xP) C¢(31) 
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tee recoswee— Sin Al sinAtya = yr? 
+ 2 [cos ¢ - sinA] sin ACxF - xa 


+ 2 [cos ¢ - sinA] cosACya - yP? 


+ 2a ’Clcos ¢ - sinA] 1/ { “tcos a= Sind ] 


2 xp sinA (xp -xa)? + 2 & Sin ACxE - Xa? C(32) 


+ 2 xF cos AClya -yFD + 2 & cosACya - yF? 


+ 2xpa’ +2 6 a’}/ Xd’ feos ein 


then our quadratic equation is, 


o= Q- 


[cos € - sin _] met ip) tees © sin 1s CC33) 


(1 - €7) +c )2¢*tcos € -sin 3 (34) 


and is of the fora, 


so that 


2 
OR =i 1 at!) 2 BCI €* >) + C CC35) 
(l= (29 9="— -2b/72 + VC 2B) - 40 
2 
2 
es fe Cif eee CD C(36) 


The following possibilities occur: 


i. 


If the term under the square root is negative then a 
solution does not exist. 
If the term under the square root is positive, and 


the value of (1 - €’) is between O and 1 then a 


solution exists; the solution is 


Chima) D2 -=B + \cB =e) 
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3. If the value of (1 - €1) obtained in step (2) is 


outside the range of O to 1 then the solution is, 


CT a Ss - \eee 


4. If the value of (1 -¢1) obtained in step (3) is 


cutside the range of 0 to 1 then no solution exists. 


Once ¢€’ has been determined and all other values have 


been found, the design is complete and the blades can be 


drawn. 


Referring to figure C2, to draw the radius on the 


leading edge the points 


x xEL *+ §&/2 cos 98 
y = yL 7- 6/2 sin 986 


are plotted over the interval of 9 between O and 180 


degrees. The radius on the trailing edge is given by, 


x xT - §&/2 cos 8 


y yT + 8/2 sin 8 


plotted over the interval of 98 between #4 and 4+180. 


Cis 


C¢(38) 


C(39) 


C640) 


To draw successive blades on the rotor the procedure is 


to rotate the axes after each blade is drawn. The procedure 


is described by equation B(SS) and BC(S6). 
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C.2 COMPUTER PROGRAM “AWEDG" 


Cee se Description 


The wedge arc blade design program, “AWEDG" 
calculates and draws wedge arc blading for an initially 
radial-outflow rotor. The difference between “AWEDG” and 
“ACIRC’ is not only the specified input parameters; "“"AWEDG" 
calculates the straight section needed on the suction side 
to constrict the flow to a know dimension at the trailing 
edge tip of the adjacent blade. The program, for simplifi- 
cation, is divided into seven parts and each part will be 
described. 

The first part (lines 1 to 460) contains the 
values of all specified parameters. On running, it also 
allows the user to specify the angles % and # to determine 
early in the program what the value of & will be. If the 
user wishes to change the value of «© he can do so. Chord 
length, camber radius and & are calculated in this section. 

The second part (lines 470 to 640) prints some 
of the output parameters and calculates the center 
coordinates of the camber line. 

In order to calculate the outlet passage size 
the preceding passage must be analyzed. The program does 
this in part three (line 650 to 800) by calculating the 


position of the trailing edge tip of the preceding pressure 
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line. Point C’ on figure C3 is calculated by rotation of 
axes. 

In part four (Clines 810 to 1020) the position of 
points C, F, T, E, and L, shown in figure C3, are 
calculated. 

Part five (beginning on line 1030) starts the 
plotting routine. The number of the blades requested to be 
drawn are specified. The pressure side radius is calculated 
and is drawn. Calculations for the straight line segments 
are performed (lines 1340 to 1570). Based on the determined 
value of the straight segment the coordinates of point I are 
dervied (lines 1630 to 2090). 

In part six the radius of the suction side is 
drawn. Plotting for the leading and trailing edge radius 
are also accomplished in part six. 

Part seven (2100 to 2140) concludes the 
procedure by printing the value of the throat diameter, 
length of the straight line segmnent and radii of pressure 


and suction circle. 


Ce2a2 Required Inputs 
The following inputs are required: 
Outer radius 
Inner radius 


Total number of blades 


angle - 
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angle % (must be consistent with ¥?-see 
discussion) | 

Diameter of the circle at the passage exit 
Leading edge and Trailing Edge thickness 
Leading edge and Trailing Edge wedge angle 


Number of blades requested to be drawn 


Listed Outputs 


The following quantities are output: 
Outer radius 

Inner radius 

Number of blades 

Angle & 

Angle @ 

Angle « 

Radial solidity 

Chord length 

Thickness to chord ratio 

Leading edge and trailing edge thickness 
Leading edge and trailing edge wedge angle 
Throat width at exit 

Radius of the pressure surface 

Radius of the suction surface 

Length of the straight section on the 


Suction side. 
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C.2.4 Program Listing 

Listed here is the program “AWEDG”™ for the 
design of wedge shaped are blades. Table C2 identifies the 
symbols used in the program with the defined parameters. 
Figure C6 and C7 illustrates computer generated shapes for 
an outer radius of ¢ inches, an inner radius of 3.2 inches 
and 60 blades. The value of ¥ and are specified 
parameters to limit the efflux angle « to between 8 and 93 


degrees. 
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TABLE Cl. 


*“AWEDG” 


AO 


A4 


AS 


BS 


BS 


CO 


Cl 


D1 


Ed 


EZ 


GCI] 


HCT] 


I1 


2 


I3 


TEXT 


List of Computer and Text Symbols 
for Wedge-Arc Blade. 


IDENTIFICATION 


efflux angle 
throat dimension (radius of 


outlet circle) 
throat diameter (diameter of 


exit circle) 
length of straight section 


length of straight section to 
be tangent to outlet circle 


c/Ro dimensionless, 
chord length 


chord length 

constant angle 
constant value 
constant value 


array for current pressure 
side x coordinates 


array for current pressure 
Side y coordinates 


counter 
Interval variable 
Interval variable 


Interval variable 


TABLE Cl. 


*AWEDGE” 


L2 


LES 


L4 


Lo 


L6G 


L7 


L8& 


MO 


M1 


M2 


M3 


M4 


MS 


List of Computer and Text Symbols for 
Wedge-Arc Blade (con’t). 


EAT 


Ye 


VE 


Yt 


ye 


¥I 


Yee 


YH 


Xe 


x 


Xt 


FE 


IDENTIFICATION 


final value on present 
pressure line, y coordinate 


final value on present 
pressure line plus blade 
thickness, y coordinate 


center point of trailing 
edge, y coordinate 


Leading edge initial suction 
side, y coordinate 


point of tangency for 
suction side, y coordinate 


point of tangency for 
suction side, y coordinate 


initial point for suction 
Side curve, y coordinate 


Leading and trailing edge 
thickness 


Leading and trailing edge 
wedge angle 


final value on present 
pressure line, x coordinate 


final value on present 
pressure line plus blade 
thickness, x coordinate 


center point of trailing 
edge, x coordinate 


Leading edge initial suction 
side, x coordinate 


TABLE Cl. 


*AWEDGE"” 


M6 


M7 


M3 


NO 


N9 


Oil 


O03 


04 


P1 


Qi 


Q2 


Q3 


Q4 


Q5 


RO 


List of Computer and Text Symbols for 
Wedge-Arc Blades (con’t). 


TEXT 


xL 


xI 


*H 


Ro 


IDENTIFICATION 


center point of leading edge 


point of tangency for suction 
side, x coordinate 


initial point for suction 
Side curve, x coordinate 


Number of blades to draw 


Number of blades to draw for 
loop 


Number of points for interval 
of arc 


Equals 1 to change # and @, 
O for no 


output print of Re 

output print of Roa 

Nt 

constant 

constant 

constant, value under square 
root term for solving for 
€1 

same as Q3 


same as Q3 


outer radius 


TABLE Ci. 


*AWEDGE” 


R1 


R2 


R3 


R4 


R6 


SO 


>2 


SZ 


S53 


SS 


TO 


T4 


T9 


X2 


X3 


X4 


X9 


List of Computer and Text Symbols for 
Wedge-Arc Blades (con’t). 


TEXT 


Ri 
e 


Re/Ro 


Rp/Ro 


Cr 


&S/Ro 


Ps 


*dPp 


Xd 


IDENTIFICATION 


inner radius 

radial chord, (Ro - Rj)/Ro 
camber line radius, 
dimensionless 

inner radius, dimensionless 
pressure radius,dimensionless 
peripheral spacing, 2nRo/Z 
radial solidity 
dimensionless peripheral 
spacing 

blade curvature with respect 
to center of rotor 

spacing angle, 3600/2 

camber angle 

loop parameter to draw arcs 
value of 1-€’ 

pressure surface center, x 
coordinate 

camber line center, x 


coordinate 


variable for x values in 
continuous loop 


variable for x values in 
rotational loop 


TABLE Cl. 


“AWEDGE” 


Y2 


Y3 


Y4 


‘3 


ZO 


List of Computer and Text Symbols for 
Wedge-Arc Blades (con’t). 


TEXT IDENTIFICATION 
Yap pressure center y coordinate 
Yd camber center y coordinate 


~ variable for y values in 
continuous loop 


= variables for y values in 
rotational loop 


Z number of blades 
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Figure C1. Relaxed Double Circular Arc Passage 
with Pressure Line Removed. 
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Figure C2. Wedge-Arc Shaped Blading Geometry. 
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Figure C3. Enlarged View of Trailing Edge. 
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Note: C’D’ia not normal to the pressure surface. 


Figure C4. Geometrical Relationship Between Throat 
Dimension and Suction Surface Radius 


Ro = 4 Ri = 3.2 No. of Blades = 60 \ 
“% = 9.17 9 = 90.68 

ehord = 0.67672 

alpha = 8.49 

radial solicvey 7. JO 3S> e727 

L.E. & T.wE. eihveknesar = 70762 

L.E. & TcE. wedge angle. —-s- 

As =» O,O0S225 

Lo = 90.0427 3506 01> 

radius of pressure circle 

radius of suction circle 


0.629980625 
0.4768272358 \) 


Figure C5. Computer Generated Shapes for a Cascade of Wedge 


Arc Shaped Baldes with 60 Blades. 
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